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ABSTRACT 

This paper examines the 48-h life cycle of a winter anticyclone occurring over North America from 18 to 20 
January 1979 using Goddard Laboratory for Atmospheres FGGE level Illb (SOP 1 ) global analyses on a 4° 
latitude by 5° longitude grid. Applying the relatively new methodology of the Zwack-Okossi equation, results 
show that anticyclonic vorticity advection and cold-air advection acted to develop the anticyclone while adiabatic 
warming in the descending air opposed development. Other forcing processes made only small contributions to 
anticyclone changes. Vertical profiles of the development quantities reveal that vorticity and temperature ad- 
vections, as well as the adiabatic wanning, maximized in the 200-300-mb layer. 


1. Introduction 

As noted recently by Tan and Curry (1993), mi- 
grating midlatitude anticyclones have been studied 
much less than have their low pressure counterparts, 
extratropical cyclones, probably due to their associa- 
tion with “fair weather” conditions. However, under- 
standing the development and movement of such an- 
ticyclones is important because of their contributions 
to the general circulation and to potential air pollution 
accumulation. Migrating anticyclones, when coupled 
with cyclones, form the transient waves of the general 
circulation, which in turn are largely responsible for 
day-to-day weather changes experienced in the midlat- 
itudes. These high pressure areas involve downward 
vertical air movement, which results in stable air and 
the absence of precipitation. Thus, when such phenom- 
ena become slow moving, or even stationary, for sev- 
eral days and occur over areas of significant pollution 
emission, air pollution problems can result. 

Wexler (1951) summarized the early meteorological 
literature pertaining to anticyclones, which he then 
classified into groups identified as cold, warm, or 
mixed. Significant radiative forcing appears to be 
largely confined to cold anticyclone development (e.g., 
Curry 1983, 1987). However, other forcing mecha- 
nisms have been noted in all anticyclone types. Among 
these are upper-level velocity convergence, identified 


* Current affiliation: Texas Natural Resources Conservation Com- 
mission, Austin, Texas. 


as a forcing mechanism itself or in association with 
anticyclonic vorticity advection (Bodurtha 1952; 
Bowling et al. 1968; Dallavalle and Bosart 1975; Boyle 
and Bosart 1983; Tan and Curry 1993) and cold-air 
advection (Fleagle 1947; Bodurtha 1952; Boyle and 
Bosart 1983; Tan and Curry 1993). 

The principal objective of this study is to add to the 
limited body of literature that currently exists regarding 
the development and decay of midlatitude migrating 
anticyclones. This objective is accomplished by stud- 
ying the synoptic-scale forcing mechanisms present 
during the life cycle of an anticyclone system that oc- 
curred over North America from 0000 UTC 1 8 January 
to 0000 UTC 20 January 1979 using a relatively new 
methodology known as the Zwack-Okossi equation 
(Zwack and Okossi 1986; Lupo et al. 1992). Indeed, 
a test of the utility of this diagnostic approach for an 
anticyclone system can be regarded as a secondary ob- 
jective of the study. This event was chosen because it 
afforded the opportunity to study a significant winter 
anticyclone that experienced both growth and decay 
over a relatively short time period. 

2. Synoptic discussion 

The horizontal and vertical structure of the anticy- 
clone system is summarized using sea level pressure, 
850-mb height and temperature, 500-mb height and rel- 
ative vorticity, and 300-mb height and wind speed 
maps. 

The anticyclone of interest here began at 0000 UTC 
1 8 January over Ontario, Canada, west of James Bay 
(Fig. la) with a central sea level pressure of 1028 mb. 
The anticyclone continued to develop and by 1 2 h later 
had moved southward across Ontario southwest of 
James Bay with a central pressure of 1037 mb (Fig. 
2a). By 0000 UTC 19 January, the anticyclone had 
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Fig 1 Observed (a) sea level pressure (4 mb, solid), (b) 850-mb heights (30 m, solid) and temperatures (5 K, dashed) 
(c, 5 TO-mb heiXs (60 m, solid) and relative vorticities (2 X 10 * s ', dashed), and (d) 300-mb he.ghts (120 m, solid) 
and wind speedL (10 ms 1 , dashed) for 0000 UTC 18 January 1979. The cross denotes the grid point with maximum 

sea level pressure. 


matured to a central pressure of 1038 mb while moving 
eastward due south of James Bay (Fig. 3a). Over the 
next 12 h, the decaying anticyclone, at 1035 mb, con- 
tinued its eastward propagation into western Quebec 
(Fig. 4a) . By the final map time, 0000 UTC 20 January, 
the anticyclone had propagated northeastward to west- 
ern Quebec just east of James Bay with a central pres- 
sure of 1029 mb (Fig. 5a). 

The surface anticyclone is associated with a double- 
ridge system in the upper-air fields, a northern ridge 
extending northeastward across the anticyclone center, 
and a southern ridge over the central and eventually 
eastern United States. At 850 mb, the temperature wave 
(Fig. lb) was out of phase with the height field, sug- 
gesting a baroclinic wave with warm-air advection west 
and cold-air advection east of the ridge system. The 
southern ridge intensified over the next 12 h (Fig. 2b) 


as it moved eastward to a position 5° west of the surface 
ridge. The thermal pattern became more out of phase 
with the height field, thus enhancing the low-level tem- 
perature advections. By 0000 UTC 19 January (Fig. 
3b), the ridge amplitude maximized over the central 
midwestem and southern United States, while the ther- 
mal pattern remained similar to the previous map time. 
During the last two map times ( Figs. 4b-5b ) , the ridge 
over the eastern United States damped as it moved to 
a position above the surface ridge. The thermal wave 
also weakened during this period. The 850-mb ridge 
height and thermal structure was similar to that reported 
by Boyle and Bosart ( 1983 ) . 

The 500-mb ridge system at 0000 UTC 18 January 
( Fig. 1 c ) lagged the surface anticyclone. During the next 
24 h (Figs. 2c -3c), the amplifying ridge moved east- 
ward to a position over the central United States, re- 
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Fig. 2. As in Fig. 1 for 1200 UTC 18 January 1979. 


maining west of the surface system. Then, over the last sart ( 1983 ) noted as a characteristic of a cold-core an- 
24 hours (Figs. 4c-5c) the ridge system damped as it ticyclone, and upstream from the jet maximum. Its po- 
moved to a position nearly above the surface anticyclone sition relative to the jet maximum, its cold core, and its 
center and ridge. Through the first four map times the relatively low tropopause (—300 mb) allow this event 
surface anticyclone features were positioned down- to be classified as a cold anticyclone, 
stream from a 500-mb relative vorticity minimum of 

about -4 X 1(T 5 s 1 . By 1200 UTC 19 January the 3 * The Zwack-Okossi equation 
vorticity minimum had shifted south of Lake Erie, and The primary diagnostic equation used in this paper is 
the anticyclone center was located just downstream from a relatively new expression known as the Zwack-Okossi 
a weaker vorticity minimum axis. The relative positions (Z-O) equation. This equation is a generalized form of 
of the 500-mb vorticity minima and the surface anticy- the Petterssen- Sutcliffe equation (Petterssen 1956, p. 
clone system were similar to those reported by Boyle 324) that diagnoses the geostrophic vorticity tendency at 
and Bosart ( 1983) and Tan and Curry ( 1993). a near-surface pressure level as forced by a set of pro- 

The 300-mb ridge locations were similar to those at cesses integrated through the entire vertical air column. 
500 mb and also followed the same amplification and Furthermore, the geostrophic vorticity tendency fields so 
decay sequence (Figs. l-5d). The other interesting determined can be relaxed to obtain near-surface geopo- 
feature to note is the location of the surface features tential height tendencies. Originally developed in quasi- 
relative to the 300-mb jet. Throughout the study period, geostrophic form by Zwack and Okossi ( 1986 ), the Z- 
the surface anticyclone center was positioned on the O equation was recently developed in both generalized 
cyclonic side of the 300-mb jet, which Boyle and Bo- and “extended 1 ’ forms by Lupo et al. ( 1992 ). 
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Fig. 3. As in Fig. 1 for 0000 UTC 19 January 1979. 
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Briefly, the development proceeds by taking the two- 
dimensional Laplacian on an isobaric surface (V 2 ) and 
Eulerian time derivative (did t) of the hydrostatic bal- 
ance equation to yield 

^^_RF‘d_VTdp 

dt dt f J p dt p 

where ^ is the geostrophic relative vorticity (g/ _1 V 2 z), 
T is the temperature, p is the pressure, R is the dry-air 
gas constant, and / is the Coriolis parameter. The sub- 
script / refers to the near-surface pressure level. Repre- 
senting d^g/dt as the difference between the actual vor- 
ticity tendency (dtydt) and the ageostrophic vorticity 
tendency {dt^ldt ) , substituting for d^ldt and dTIdt 
using the vorticity equation and first law of thermody- 
namics, respectively, and integrating from p t to an upper- 
level p, that is sufficiently high to include most of the 
atmospheric mass while still providing usable data yields 
the generalized expression [see Eq. (9) in Lupo et al. 
1992] . As noted in Lupo et al., it is appropriate for stud- 


ies of large-scale circulations to simplify the generalized 
expression based on scale analysis and many sample cal- 
culations of individual terms. This yields the “ex- 
tended” form of the Z-O equation 



(B) (C) (D) 


+ Pd 


f 


k* V x F dp , 
(E) 


(2) 


where Pd = (pi - p,)”\ Ci absolute vorticity, V is 
the horizontal wind vector, u is the vertical motion in 
isobaric coordinates ( dpfdt ), Q is the diabatic heating 
rate, F is the frictional force, and S is the static stability 
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Fig. 4. As in Fig. 1 for 1 200 UTC 19 January 1979. 


parameter [= -(T/9)(d0/dp)], 0 is potential temper- 
ature). Equation (2) is so-named following the termi- 
nology of Tsou et al. ( 1987) in their treatment of the 
height tendency equation. The implication is that the 
quasigeostrophic form is extended to include the influ- 
ence of nonquasigeostrophic forcing in the right-hand- 
side terms. In all terms of (2) the integrand of the p } - 
p , integral represents the horizontal divergence forced 
by each process at level p. Term A (VADV) represents 
local divergence (convergence) that is produced as the 
atmosphere readjusts to a spatially varying, geostroph- 
ically balanced state in response to local cyclonic ( an- 
ticyclonic) vorticity advection. Terms B (TADV), C, 
and D (ADIA) represent heating (cooling) processes 
that force horizontally nonuniform upward (down- 
ward) displacements of the pressure surfaces above the 
heating (cooling) maximum, resulting in local hori- 
zontal divergence (convergence). Term E (ERIC), 
which is restricted to the boundary layer in this paper, 
represents the divergence or convergence forced by 


frictional processes. Term C includes the effects of la- 
tent heat release (LATH), infrared radiation (IRCH), 
and boundary layer sensible heating ( SENS ) . Note that 
the Z-0 equation explicitly includes both dynamic and 
thermal forcing processes at all levels above the de- 
veloping surface feature. 

4. Computational methods 

a. Data and finite differencing 

The Goddard Laboratory for Atmospheres (GLA) 
6-h interval special observation period 1 level Illb 
analyses of the FGGE level lib dataset (Baker 1983) 
were used for this research. This study utilized as its 
computational domain the portion of the original 
GLA global analyses extending from 22° to 82°N and 
0° to 120°W. This domain covers approximately all 
of North America and adjacent portions of the North 
Atlantic Ocean. Mandatory-level data were interpo- 
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Fki. 5. As in Fig. 1 for 0000 UTC 20 January 1970. 


lated or extrapolated to standard isobaric surfaces strophic vorticity tendencies with zero initial and lateral 
from 1050 to 50 mb in 50-mb increments. The hor- boundary values. Relaxation was terminated when the 

izontal resolution is 4° latitude by 5° longitude. absolute difference between two successive iterations 

The results presented are confined to moving do- was less than 10 ' m s 1 at all grid points for the 
mains extending over 32° latitude and 35° longitude VADV, TADV, ADIA, and LATH terms and 10 
in order to focus on the evolution of the surface an- ms 1 for the remaining terms. Relaxation of the ther- 
ticyclone previously described. The vertical compu- modynamic terms followed the application of the hor- 
tational domain extends from the near-surface level izontal smoothing operator discussed in Lupo et al. 
(p ; ), defined as the first 50-mb layer above the sur- ( 1992). In order to identify the tendency quantities 
face (always either 1000 or 950 mb), to an upper with development (i.e., where the propagation contri- 
limit of 100 mb. The near-surface level was so de- bution to the tendency would be expected to be small), 
fined in order to specify a pressure level that was results are also summarized for the anticyclone center, 
always above the surface but was sufficiently close Vertical motions were determined using an extended 
to the surface to be representative of surface anticy- form of the omega equation (see Tsou et al. 1987; Lupo 
clogenesis. et al. 1992) with u> set to zero at the surface and 50 mb. 

The Z-O equation calculations applied second-order 
finite differencing to determine all derivatives and the Diabatic heating and frictional processes 
trapezoidal rule for the evaluation of vertical integrals. 

Geopotential height tendencies were obtained by two- Latent heat release, sensible heating, and boundary 
dimensional sequential overrelaxation of the geo- layer frictional processes were obtained using the pro- 
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cedures described in Lupo et al. (1992). Infrared ra- 
diation was determined using Sasamori's (1968) 
scheme with a water vapor emissivity correction sug- 
gested by Ramanathan et al. (1983). This scheme, 
which contains the effects of water vapor, carbon di- 
oxide, and ozone, determines the longwave cooling 
rates as the difference between integrated upward and 
downward radiative fluxes in a cloudless atmosphere. 
To account for clouds, a randomly overlapping scheme 
developed by Harshvardhan et al. ( 1987) was utilized. 
A cloud was assumed to exist at any point that had a 
latent heat value greater than zero. 

c. Comparison procedures 

Mean absolute values ( MA V ) and correlation coef- 
ficients (CC) were employed to statistically compare 
calculated Z-O with observed geopotential height ten- 
dencies. The observed tendencies were estimated as 1 2-h 
finite differences centered on the period from 6 h before 
to 6 h after each map time. The mean absolute value is 
the mean of the absolute value of the height tendencies 
over all grid points. All statistics were computed for 
the anticyclone domains represented in Figs. 6 and 7. 
These statistics were also averaged over the complete 
time period from 0000 UTC 1 8 January to 0000 UTC 
20 January 1979. 

5. Comparison of observed and Z-O tendencies 

An integrated view of the comparability of ob- 
served near-surface geopotential height changes 
(OBS) with Z-O equation values (Z-O) using 
MAV and CC statistics is given in Table 1. Results 
indicate that the observed tendencies averaged over 
the 48-h period are approximately 17% larger in 
magnitude than the Z-O tendencies at the near-sur- 
face level. The larger percentage differences at 0000 
UTC 19 January occurred because of underestimates 
of Z-O height falls in the wake of the anticyclone in 
the western portion of the domain. At 0000 UTC 20 
January, the difference can be largely attributed to 
the placement of the Z-O height rise maximum one 
grid interval too far east and thus partially outside of 
the domain represented in the Table 1 statistics. In 
neither case is the anticyclone diagnosis compro- 
mised. Furthermore, the correlation coefficients 
(5*0.786) suggest good pattern comparability be- 
tween the two fields at all times. 

To examine both the magnitude and pattern com- 
parability further, comparisons of the two fields are 
made for 1200 UTC on both 18 and 19 January (Fig. 
6). These map times were selected because they cor- 
respond to the midpoint of both the development and 
decay phases of the surface anticyclone. At 1 200 UTC 
18 January, the near-surface OBS and Z-O patterns 
and magnitudes are very similar, with a particularly 
good match in the area of the height increases over the 


eastern Great Lakes. The near-surface OBS and Z-O 
tendencies patterns are also quite similar at 1200 UTC 
19 January. The OBS and Z-O maxima are located 
over the Atlantic Ocean at about 40°N, 65°W with an 
OBS maximum that is 28% larger than the Z-O. The 
patterns and magnitudes are sufficiently similar to con- 
clude that the Z-O model provides reliable height ten- 
dency estimates for the surface anticyclone. 

6. Diagnostic results 

a. Basic forcing processes 

In the interest of brevity, the discussion in this sec- 
tion is restricted to the two map times noted above. 
Furthermore, maps of individual forcing contributions 
(Fig. 7) are restricted to those that most consistently 
dominate the anticyclone evolution, the horizontal vor- 
ticity and temperature advections, and the adiabatic 
warming. Other forcing processes are examined at the 
anticyclone center in the next two sections. 

At 1200 UTC 18 January, midway through the pe- 
riod of anticyclone pressure rise, an area of height in- 
creases encompassed the surface center, which was 
northwest of the total tendency maximum (see upper- 
right panel in Fig. 6). Figure 7 (left-hand panels) re- 
veals that these height increases were largely forced by 
integrated anticyclonic vorticity and cold-air advec- 
tions. The combination of these two processes was suf- 
ficient to offset the influence of adiabatic warming 
found within the descending air associated with the an 
ticyclone. During the anticyclone decay period, as rep- 
resented by the 1200 UTC 19 January field ( lower right 
panel in Fig. 6), the surface anticyclone center was 
positioned near the zero tendency isopleth, again north- 
west but farther removed from the total height rise max- 
imum. 

b. Anticyclone center averages 

Recall that anticyclone center values (see section 4a) 
are used to isolate the development component of the 
height tendency. In order to summarize the relative 
contribution of the different terms, a quantity referred 
to as the “development sum,” formed by adding the 
terms that contribute to anticyclone development ( pos- 
itive height tendencies), was computed for the 48-h 
average statistics. All terms, including those that re- 
sulted in negative tendencies, were then evaluated rel- 
ative to this sum. 

The 48-h average total height tendency (ZTOT) for 
the anticyclone center (Fig. 8, top) was 0.69 X 1 0 ' 3 
m s “ 1 . The development was clearly dominated by hor- 
izontal anticyclonic vorticity and cold-air advections, 
which contributed 68% and 27%, respectively, to the 
development sum of 2.37 X 10~ 3 m s 1 . The primary 
balancing process was adiabatic warming, which was 
59% of the development sum. All other terms were 
small, reflecting the small values typical of these quan- 



2280 


MONTHLY WEATHER REVIEW 


Volume 123 



Fig. 6. Near-surface observed (left) and Z-O (right) height tendencies for 1200 UTC 18 January 
and 1200 UTC 19 January 1979. Increment and units are 0.50 x 10 'ms 1 . 


tides at the anticyclone center and/or the small Lapla- 
cians to be expected from their nearly uniform hori- 
zontal distributions. 

Examining the two individual map times previously 
discussed, the midpoint of vigorous anticyclone de- 
velopment (1200 UTC 18 January) was marked by 
large height increases (2.21 X 10 3 m s 1 ). The most 
significant forcing terms ( VADV, TADV, and ADI A) 
were also large at this time. At 1200 UTC 19 January 
a small positive total tendency was evident (0.54 
X 10" 3 m s -1 ), due primarily to closer balance be- 
tween reduced adiabatic warming and anticyclonic 
vorticity and cold-air advections. Failure to obtain a 
negative value for the total tendency at this time is 
most likely due to errors in the calculations of the 
individual forcing terms, which are more acutely felt 
when the grid point of interest is near the zero ten- 
dency isopleth. 


c. Vertical profiles 

Clearly, much of the forcing of near-surface height 
tendencies occurs as a result of processes well above 
the surface. These contributions are examined for the 
two map times previously discussed using vertical pro- 
files of quantities contributing to each forcing term at 
the anticyclone center. 

Figure 9 shows vertical profiles of relative vorticity 
and absolute vorticity advection. The relative vorticity 
during anticyclone development ( 1 8 January ) was pos- 
itive above 750 mb with a maximum at 400 mb. During 
the decay phase ( 19 January ), the positive relative vor- 
ticity values were smaller with the maximum elevated 
to 200 mb. Negative values extended through a deeper 
layer from the surface to 500 mb. The cyclonic vortic- 
ities evident in the upper portion of these profiles can 
be explained by recalling the position of the anticy- 
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Fig. 7. Contribution to the total Z-O height tendency due to horizontal vorticity advection 
(upper), horizontal temperature advection (middle), and adiabatic heating /cooling (lower) for 1200 
UTC 18 January (left-hand panels) and 19 January 1979 (right-hand panels). Increment and units 
are 0.50 X 10 3 m s - ’. 


clone relative to the upper-air flow. Examining Figs, tantly, throughout the period was located on the cy- 
1 -5d, one sees that the anticyclone was initially located clonic side of the jet. The vertical profiles of vorticity 
under a cyclonic curvature region and, more impor- advection show anticyclonic vorticity advection 
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Table 1. Near-surface mean absolute values (MAV, 10 'ms ') 
and correlation coefficients (COR) for 0000 UTC 18 January 
through 0000 UTC 20 January. AVG represents the value aver- 
aged over the five map times. DIFF represents the percentage 
difference between OBS and Z-0 MAV values, that is, 
100{ lOBSMAV - (Z-OMAV))/OBSMAV }. 


MAV 



N 

1 

O 

OBS 

DIFF 

COR 

0000 UTC 18 January 

0.747 

0.783 

5 

0.943 

1200 UTC 18 January 

0.895 

0.957 

6 

0.935 

0000 UTC 19 January 

0.773 

1 .042 

26 

0.786 

1200 UTC 19 January 

0.880 

0.938 

6 

0.892 

0000 UTC 20 January 

0.705 

1.081 

35 

0.91 1 

48-h AVG 

0.800 

0.960 

17 

0.893 


(AVA) throughout the column during development, 
with a maximum at 300 mb, similar to the results of 
Dallavalle and Bosart ( 1975). By 19 January the AVA 
had decreased markedly, with weak positive values ap- 
parent above 250 mb. The temporal changes in AVA 
that were evident over the anticyclonic center are con- 
sistent with a tilted baroclinic wave. As the wave de- 
veloped and then decayed, the upper-air ridge and the 
area of minimum vorticity in the upper troposphere 
overtook the surface anticyclone, resulting in a de- 



Fig. 8. Total and individual term contributions for the Z-O equa- 
tion height tendency (10“ 3 m s ') at the anticyclone center for the 
48-h average (top), 1200 UTC 18 January (middle), and 1200 UTC 
19 January 1979 (bottom). 



Fig. 9. Vertical profiles of relative vorticity (10 5 s solid) and 
horizontal vorticity advection (10“ 9 s' 2 , dashed) for 1200 UTC 18 
January and 1200 UTC 19 January 1979 at the anticyclone center. 
Ordinate is pressure in millibars. 


crease in vertically integrated AVA above the anticy- 
clone. 

Figure 10 shows vertical profiles of temperature ad- 
vection for the anticyclone center. Similar to the vor- 
ticity advection, the temperature advection profiles 
show the forcing due to cold-air advection to be strong- 
est in the upper troposphere /lower stratosphere. These 
results are similar to those of Bodurtha (1952) and 
Fleagle ( 1947). The profiles also show during anticy- 
clogenesis (1200 UTC 18 January) a secondary cold- 
air advection maximum in the boundary layer. By 1 200 
UTC 19 January, decreased upper-level cold-air advec- 
tion was accompanied by lower-troposphere warm-air 
advection, thus reducing the integrated cold-air advec- 
tion. 

Figure 1 1 displays vertical profiles of the product of 
the static stability and vertical motion (Su;). Note that 
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Fig. 10. As in Fig. 9 with vertical profile of horizontal 
temperature advection (10 4 K s ’). 


these profiles are opposite in sign to the profiles for the 
temperature advection. During anticyclogenesis Suj ex- 
hibits a maximum at 300 mb, where increasing S values 
[S increases rapidly above 300 mb, see Lupo et al. 
( 1992, Fig. 23)1 are coupled with sufficiently strong 
downward motions. During the decay phase, weakened 
vertical motions resulted in smaller Su) values. These 
profiles indicate subsidence and adiabatic warming 
over the anticyclone that decreased near-surface height 
tendencies throughout the development and decay of 
the anticyclone. 

The vertical profiles of friction and sensible heating 
(not shown) yield nonzero values below 850 mb only, 
with friction and sensible heating yielding tendencies 
in accordance with the results shown in Fig. 8. Profiles 
of latent heat release and infrared cooling (also not 
shown) are zero for the former and are consistently 
negative for the latter, but they are an order of mag- 
nitude smaller than the temperature advection through- 
out the column. 


7. Discussion 

A diagnosis of a winter North American anticyclone 
that occurred from 0000 UTC 1 8 January to 0000 UTC 
20 January 1979 was performed using the near-surface 
height tendencies derived from the Zwack-Okossi 
equation. 

Using horizontal fields and vertical profiles, the re- 
sults demonstrate that anticyclonic vorticity advection, 
dominated by upper-tropospheric- lower-stratospheric 
advections, was the most significant contributor to the 
total near-surface height increases found upstream, 
within, and downstream from the anticyclone center. 
Temperature advection, dominated by cold-air advec- 
tion at the same upper levels, was the next most sig- 
nificant development term. Thus, dynamic and ther- 
modynamic transport processes in the upper tropo- 
sphere and lower stratosphere had a profound influence 
on the anticyclone’s development. This is consistent 
with the results of Austin and Shapiro ( 1951 ) and Bo- 



Fig. 1 1. As in Fig. 9 with vertical profile of adiabatic 
warming (Su, 10 4 K s '). 
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durtha ( 1975), who found that mass increases associ- 
ated with anticyclone pressure increases are most 
strongly contributed by the upper troposphere/lower 
stratosphere. Further, the dominance of the two advec- 
tion processes for the modest anticyclone studied here 
corresponds to similar results to be inferred from Tan 
and Curry ( 1993 ) for a very strong winter anticyclone. 
The adiabatic term, because of adiabatic warming in 
the descending air, contributed negatively to the de- 
velopment of the anticyclone and was the same order 
of magnitude as the vorticity and temperature advection 
terms. Boundary layer sensible heating, latent heat re- 
lease, boundary layer friction, and infrared cooling 
made only minor contributions to anticyclone devel- 
opment. The latter, perhaps surprising, result is simi- 
lar to that found by Tan and Curry (1993), both of 
which stand in contrast to the model results of Curry 
(1987). 

Finally, comparison with the results of Lupo et al. 
( 1 992 ) reveals that aside from the absence of a latent 
heat release influence, the large-scale forcing mecha- 
nisms important in the development of this migrating 
anticyclone case are the same, although of opposite 
sign, as those for migrating cyclones. Horizontal vor- 
ticity and temperature advection, both of which maxi- 
mize at upper levels, and adiabatic heating (cooling) 
are dominant forcing processes for both phenomena. 
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